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利用 端点 延 拓 提高 LS+NN 模型 的 UT1-UTC 预报 精度 
RAT, ER 
(1. 宝鸡 文理 学 院 电子 电气 工程 学 院 ， 陕 西 宝鸡 721013 
2. 中 国 科 学 院 国家 授时 中 心 ， 陕 西 西安 710600) 

摘 E: 现 有 UT1-UTC 预 报 模式 在 进行 周期 项 与 残 差 项 拟 合 分 离 时 ， 通 常 没有 考虑 最 小 
ZR (least-squares, LS) 拟 合 序 列 的 端 部 效应 ， 预 报 精度 难以 取得 较 大 提高 。 针 对 LS 拟 
合 存 在 的 端 部 效应 ， 首 先 采 用 灰色 模型 在 UT1-UTC 序 列 的 两 端 进行 数据 延 拓 ， 形 成 一 个 新 序 
列 ， 然 后 对 新 序列 进行 LS 拟 合 ， 最 后 再 联合 LS 外 推 和 神经 网 络 (LS+NN ) 模型 对 UT1-UTC 原 始 
序列 进行 外 推 。 结 果 表 明 ， 对 UT1-UTC 序 列 进 行 端点 数据 延 拓 再 进行 LS 拟 合 ， 能 够 有 效 地 改 
善 LS 拟 合 序 列 的 端 部 效应 ; 相对 于 常规 LS+NN 模 型 ， 端 部 效应 改善 的 LS+NN 模 型 的 UT1-UTC 预 
报 精 度 有 一 定 提高 ， 尤 其 对 中 长 期 预报 精度 提高 更 为 明显 。 
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0 引言 


HESR É 
括 岁差 、 章 动 、UT1-UTC 及 极 移 的 两 个 分 量 x 与 y, EOP 是 实现 地 球 参考 坐标 系 与 天 球 参 


考 坐 标 系 相互 转换 的 必要 参数 , 在 天 文 地 球 动力 学 研究 、 卫 星 导 航 及 深 空 探测 与 跟踪 等 领域 
均 有 重要 作用 。 其 长 干涉 基线 测量 (Very Long Baseline Interferometry, VLBI) 和 全 球 卫 
星 导航 系统 (Global Navigation Satellite System, GNSS) 等 空间 测 地 技术 是 测量 EOP 的 主 
要 手段 ， 其 中 UT1-UTC 的 测量 精度 为 3~5hs， 但 复杂 的 资料 处 理 过 程 致使 EOP 的 获取 存在 
一 定 的 延迟 。 由 于 深 空 探测 与 跟踪 等 领域 对 EOP 实时 观测 值 有 重要 需求 , 而 卫星 自主 导航 又 
对 EOP 中 长 期 预报 值 有 重要 需求 , 因此 对 EOP 进行 高 精度 的 短期 与 中 长 期 预报 具有 重要 的 现 
实意 义 。 在 EOP 的 几 个 分 量 中 ，UT1-UTC 是 变化 最 快 、 最 难 精确 预报 的 一 个 分 量 。 
当前 预报 UT1-UTC 的 方法 有 多 种 , 其 中 多 数 方法 联合 最 小 二 乘 (least-squares，LS) 外 推 和 
神经 网 络 (neural network，NN) 或 自 回 归 (autoregressive，AR) 等 模型 进行 预报 *”， 这 些 
方法 的 思路 为 ， 首 先 采 用 最 小 二 乘 外 推 模型 提取 UT1-UTC 序 列 中 的 周期 项 ， 并 对 周期 项 进行 
外 推 ， 然 后 利用 神经 网 络 或 AR 等 模型 等 对 LS 拟 合 残 差 进行 建 模 、 预 报 ， 最 后 再 将 周期 项 和 残 
差 项 的 外 推 值 相 加 以 得 到 UT1-UTC 预 报 值 。 在 预报 实践 中 发 现 , 在 利用 LS 外 推 模型 对 UT1-UTC 
观测 资料 进行 拟 合 时 , 在 拟 合 序列 的 两 端 存在 发 散 畸 变现 象 , 这 种 现象 在 数据 处 理 中 称 为 端 
部 效应 ””。 端 部 效应 使 残 差 项 、 周 期 项 的 预报 值 出 现 偏差 , 最 终 导致 UT1-UTC 预 报 值 不 准确 。 
本 文 针 对 UT1-UTC 预 报 中 LS 拟 合 出 现 的 端 部 畸变 现象 ,在 采用 LS 外 推 模型 对 UT1-UTC 观 测 资 
料 进行 拟 合 之 前 , 先 利用 灰色 模型 对 UT1-UTC 观 测 资 料 进行 数据 延 拓 ， 即 在 UT1-UTC 序 列 的 两 
端 增加 应 用 灰色 模型 延 拓 的 若干 数据 点 ,形成 一 个 新 序列 ,然后 利用 LS 外 推 模型 对 新 序列 进 
行 拟 合 ， 最 后 再 联合 LS 外 推 模型 和 神经 网 络 (LS+NN) 对 原始 UT1-UTC 序 列 进行 外 推 ， 这 样 就 
可 以 将 LS 拟 合 出 现 的 端 部 畸变 现象 移 至 新 序列 的 两 端 。 数 值 分 析 表 明 ， 通 过 在 UT1-UTC 序 列 
两 端 增加 统计 延 拓 数据 ， 可 以 有 效 地 抑制 LS 拟 合 出 现 的 端 部 效应 现象 ， 并 且 可 以 显著 提高 
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转运 动 可 以 用 地 球 定 向 参数 (Earth orientation parameters，EOP) 来 表征 。EOP 包 
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UT1-UTC 预 报 精 度 。 


1 端点 延 拓 算法 

灰色 模型 可 以 在 系统 变化 规律 不 明确 情况 下 利用 少量 的 已 知 信息 来 推 知 系统 的 发 展 趋势 ， 
且 预 测 精 度 高 、 建 模 简 单 。 一 方面 ， 灰 色 模 型 可 以 利用 较 少 的 观测 数据 进行 UTI1-UTC 序 列 的 
端点 外 推 延 拓 ， 所 外 推 的 的 序列 可 以 在 一 定 程度 上 反映 UT1-UTC 序 列 的 整体 变化 趋势 ”， 另 
一 方面 ， 基 于 端点 延 拓 算法 所 形成 的 UT1-UTC 新 序列 仅 用 来 求解 LS 模型 系数 ， 目 的 是 将 LS 拟 
合 出 现 的 端 部 效应 移 至 新 序列 的 两 端 ， 而 两 端 序列 并 不 参与 UT1-UTC 预 报 ”"。 因 此 ， 基 于 灰 
色 模 型 的 端点 延 拓 算法 用 于 改善 LS 拟 合 的 端 部 效应 是 完全 可 行 的 。 下 面 曾 述 基于 灰色 模型 的 
端点 延 拓 算法 。 

灰色 模型 建 模 的 基本 思路 为 ， 首 先 利 用 少量 原始 观测 值 通过 累加 生成 操作 (accumulated 
generating operation, AGO) 生成 一 次 累加 序列 ， 然 后 利用 此 累加 序列 建立 相应 的 离散 预 
测 模型 ， 并 获得 累加 序列 的 预测 值 ， 最 后 再 通过 累 减 生成 操作 (inverse accumulated 
generating operation, IAGO) 获得 原始 序列 的 预测 值 。 下 面 直接 给 出 灰色 预报 模型 GM(1, 1) 


ja 


_ 的 一 般 形式 : 

= (kt p)=[x(1l) ulaj “? 1 e") (1) 
(O 式 中 ，U 、2Q 为 灰色 模型 系数 ;天 为 参与 建 模 的 原始 观测 值 个 数 ， 万 为 预报 跨度 。 

人 

a BEX = {x(1), x(2),-++, x(n) } 为 原始 离散 序列 ， 则 灰色 模型 端点 延 拓 算法 如 下 : 


oo D 对 于 原始 离散 序列 天 的 左 端 ， 在 左 端 点 处 取 大 个 观测 值 ， 其 反 序列 为 
X, = {x(k), x(k-1),---, x(1)} (2) 


将 天 | 作为 左 端 建 模 序 列 进行 灰色 预测 ， 得 到 预测 序列 ; 


X, = {x(0), 2C), =, x(1— p)} (3) 


则 瑟 , 的 反 序 列 ， 即 原始 离散 序列 到 的 左 端点 延 拓 序 列 为 


X, = {x(1— p), X(2— p), +, x(0)} (4) 


2) FRR BOP XA Ain, ZEA ke SME, EFA A: 


X,= {x (n—k +1), xO(n—k +2), e, x(n) } (5) 


将 羡 , 作 为 右 端 建 模 序列 进行 灰色 预测 ， 得 到 预测 序列 ; 
X,= {x(n+1), x(n+2),.…, X(n+ p)} (6) 


将 序列 着,、 关 ;分别 作为 原始 离散 序列 兰 的 无 右 端 点 延 拓 值 ， 得 到 延 拓 后 的 离散 序列 : 


X’ = {x(1—- p),……, X(0), x(1),:…, x(n), xX(n+1),.…, X(n+ p) (7) 
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2 预报 模型 
2.1 最 小 二 乘 外 推 模型 
UT1-UTC 序 列 中 的 长 期 趋势 项 、 周 年 项 及 半年 项 利用 如 下 模型 进行 拟 合 : 


2 
f= + t+ |d cos D+esin( | (8) 
i=l 
式 中 ， ，、 ,表示 周年 项 和 半年 项 的 振荡 角 频 率 ， 取 ,=2 /365.24 、 


,=2 /182.62®; 、 、4d, 和 6 为 未 知 参数 ， 可 通过 最 小 二 乘法 求解 。 


2.2 神经 网 络 模型 
提取 长 期 趋势 项 和 周期 项 后 的 LS 拟 合 残 差 序列 中 含有 剩余 的 短 周期 成 分 , 也 含有 非 线性 成 
分 ， 对 于 这 些 成 分 ， 利 用 神经 网 络 进行 建 模 较 为 合理 。 对 于 LS 拟 合 残 差 序 列 


fe(t),t=1,2,--,n}> (ERMA (Ae 1、t 2、.…、1 mm 时 刻 的 历史 


REWE (t D E 2)、…， ( m) 来 进行 预测 ， 则 可 建立 从 输入 模式 到 输出 模式 的 
数据 映射 f R” 一 RR， 预 测 模型 可 表示 为 


e(t)= f([e(t—1),e(t—2),:…, E(t—m)]) (9) 


AP, 11 为 和 代 入 维 数 ， 本 文采 用 极限 学 习 机 (extreme learning machine, ELM) 算法 训练 神 
经 网 络 名 ， 从 而 建立 从 从 输入 模式 到 输出 模式 的 数据 非 线性 映射 。 在 残 差 预报 过 程 中 ， 当 已 
经 预报 出 一 部 分 残 差 值 的 情况 下 , 后 面 的 残 差 是 根据 其 所 对 应 的 历 元 通过 移动 窗口 的 方法 而 
获得 中。 
2. 3 预报 流程 

本 文 模型 与 常规 LS+NN 模 型 的 差别 之 处 在 于 ， 本 文 模型 在 对 UT1-UTC 序 列 进行 LS 拟 合 之 前 ， 
先 应 用 灰色 模型 在 UT1-UTC 序 列 两 端 进行 数据 延 拓 ， 然 后 再 对 延 拓 后 的 UT1-UTC 序 列 进行 LS 
拟 合 ， 建 立 周 期 项 和 趋势 项 LS 外 推 模型 ， 最 后 利用 LS+NN 模 型 对 原始 UT1-UTC 序 列 进行 预报 。 
值得 说 明 的 是 ， 在 对 UT1-UTC 序 列 进行 外 推 之 前 ， 需 要 根据 国际 地 球 自 转 与 参考 系 服务 
(International Earth Rotation and Reference Systems Service, IERS) 协议 扣除 UT1-UTC 
序列 中 周期 为 5d 一 18. 6a 的 62 个 固体 地 球 带 谐 潮 汐 项 ”， 同 时 还 必须 扣除 疼 秒 。 将 本 文 提出 
的 模型 称 之 为 端 部 效应 改善 的 LS+NN Cedge-effect corrected LS+NN, ECLS+NN) PAm, H 
预报 流程 如 图 1 所 示 。 


A E = / 序列 
图 1 ECLS+NN 模型 的 预报 流程 
Fig.1 Flowchart of the ECLS+NN prediction model 


3 模型 验证 
3.1 数据 来 源 
本 文 试 验 所 用 的 UT1-UTC 数 据 来 自 IERS 发 布 的 EOP 08 C04 序 列 ， 时 间 足 度 为 2000-1-1 一 


2016-10-1， 采 样 间隔 为 1d。 
3. 2 精度 评定 标准 
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为 了 衡量 UT1-UTC 预 报 结果 的 精度 ， 选 用 平均 绝对 误差 (mean absolute error, MAE) 作 
为 精度 评定 标准 ， 其 计算 公式 为 


1 N 
MAE, =— |F” o| (10) 
N 


i=l 


ah, NAAR FY. OP 分 别 为 ;点 的 预报 值 和 观测 值 。 


结果 分 析 


1.1 一 2016.9.20 期 间 的 UT1-UTC 观 测 资料 进行 试验 分 析 ， 其 中 2010.1.1 一 
2016. 9. 20 为 预测 期 ， 建 模 数据 长 度 为 10a， 每 隔 7d 预 报 1 次 ， 总 共 进 行 了 300 次 预测 。 
为 了 验证 端点 延 拓 方法 对 LS 拟 合 端 部 效应 的 改善 效果 ， 首 先 对 比 端 点 延 拓 前 后 LS 拟 合 效 


果 。 图 2 分 别 绘 出 了 端点 数据 延 拓 前 后 LS 拟 合 序 列 左 、 右 端 360 个 历 元 的 拟 合 残 差 数据 点 ， 拟 


合 时 段 为 2002. 1. 1 一 2011. 12. 31， 其 中 ， 端 点 数据 延 拓 个 数 为 240， 即 在 UTI-UTC 观 测序 用 列 的 
左 、 右 两 端 各 延 拓 120 个 数据 点 。 从 图 2 可 以 看 到 ， 与 直接 对 UT1-UTC 原 始 上 P 


基于 灰色 模型 端点 延 拓 后 的 UT1-UTC 序 列 在 左 、 右 两 端 LS 拟 合 的 残 差 值 Hh, 这 说 明 端 点 
拓 方 法 可 以 提高 LS 模型 在 UT1-UTC 序 列 左 、 右 两 端的 拟 合 准确 性 ， We 


的 端 部 效应 现象 。 
a4 (a) UT1-UTC 序 列 左 端 拟 合 残 差 Par (b) UT1- UTC 序列 右 端 拟 合 残 差 
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图 2 LS 和 ECLS 模型 在 UT1-UTC 序列 左 、 右 两 端的 拟 合 残 差 
Fig.2 LS and ECLS fitting residuals of the UT1-UTC time-series at boundaries 


为 了 进一步 检验 基于 灰色 模型 的 端点 延 拓 方 法 对 LS+NN 模 型 预报 UT1-UTC 的 改善 效果 , 本 文 


分 别 利 用 LS+NN 模 型 与 ECLS+ 
结果 的 MAE 曲 线 , 表 1 统 计 了 LS+NN 模 型 和 ECLS+NN 模 型 在 不 同 预报 长 
度 下 的 MAE 值 。 从 图 3 与 表 1 可 


ECLS+NN 模 型 UT1-UTC 预 报 


NN 模 型 对 UT1-UTC 进 行 1 300d 的 预报 ， 图 3 绘 出 了 LS+NN 模 型 和 


以 发 现 ，ECLS+NN 模 型 相对 于 LS+NN 模 型 的 UT1-UTC 预 报 精 度 在 各 


个 跨度 都 有 不 同 程度 的 提高 ， 其 中 ， 对 于 路 度 为 1 一 90d 的 预报 ， 精 度 提高 在 10% 以 内 ， 从 第 
90d 开 始 ，ECLS+NN 模 型 的 预报 精度 提高 越 来 越 显 车 ， 最 大 提高 达 40%， 这 说 明 ECLS+NN 模 型 对 
于 UT1-UTC 中 长 期 预报 精度 改进 更 为 显著 。 
表 1 LS+NN 模 型 和 ECLS+NN 模 型 不 同 跨度 UT1-UTC 预 报 结果 的 MAE 统 计 
Tab.1 MAE statistics of the UT1-UTC predictions by the LS+NN and ECLS+NN models 


跨度 /d 


LS+NN/ms — ECLS+NN/ms 精度 改善 百分比 /% 


1 
5 


0.088 0.081 8 
0.399 0.394 1 
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10 0.962 0.952 1 
20 2.214 2.199 1 
30 3.475 3.344 4 
60 6.427 6.190 4 
90 10.695 9.583 10 
120 17.067 12.822 25 
180 29.344 23.925 18 
240 40.915 26.617 35 
270 48.384 28.866 40 
300 54.359 33.408 39 
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图 3 LS+NN 模型 和 ECLS+NN 模型 UT1-UTC 预报 结果 的 MAE 比较 
Fig.3 MAE comparison between the LS+NN and ECLS+NN models forthe UT1-UTC predictions 

4 结果 分 析 

提出 了 一 种 利用 端点 延 拓 改善 LS 拟 合 序 列 端 部 效应 的 方法 ， 该 方法 在 对 UT1-UTC 序 列 进行 
LS 拟 合 之 前 ， 首 先 在 UT1-UTC 序 列 两 端 增加 统计 延 拓 数据 ， 然 后 再 对 数据 延 拓 后 的 新 序列 进 
行 LS 拟 合 , 目的 是 将 LS 拟 合 存在 的 端 部 畸变 移 至 新 序列 的 两 端 ， 从 而 抑制 原始 拟 合 序列 的 端 
部 畸变 。 试 验 结果 表明 ， 通 过 在 UT1-UTC 观 测序 列 的 两 端 增加 用 灰色 模型 延 拓 出 的 外 推 数据 
点 ， 然 后 再 进行 LS 拟 合 ， 能 够 有 效 地 抑制 端 部 效应 的 影响 ， 与 常规 LS+NN 预 报 模型 相 比 ， 
ECLS+NN 模 型 的 UT1-UTC 预 报 精度 在 各 种 跨度 均 有 不 同 程度 提高 , 其 中 在 跨度 为 中 长 期 预报 中 
精度 提升 尤为 明显 ， 因 此 ECLS+NN 模 型 更 适用 于 UT1-UTC 中 长 期 预报 。 
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Improving the performance of the LS+NN model for 


UT1-UTC forecast with the edge extension 
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Abstract: The prediction accuracy of UT1-UTC can be easily affected by the edge distortion 
of least-squares (LS) fitting time-series, referred to as the edge effect in data-processing fields, 
when periodic oscillations and residuals are separated by LS fitting. In order to alleviate the edge 
effect, the original UT1-UTC time-series is first extended on both boundaries using the GM(1,1) 
grey model in this work. A LS extrapolation model is then set up using the extended time-series to 
remove the edge distortion to the boundaries. Finally, UT1-UTC predictions are generated by the 
combination of the edge effect correlated LS (ECLS) extrapolation model and extreme learning 
machine neural network (ECLS+NN). The numerical experiments show that the edge effect can be 
remarkably alleviated with the presented approach. In addition, the accuracy of the UT1-UTC 
short-term predictions by the ECLS+NN method is better than that obtained by the LS+NN 
approach, but only slighter. The medium- and long-term predictions, however, are noticeably 
more accurate than those by the LS+NN solution. 

Key words: Universal time (UT1); Prediction; Edge extension; Grey model; 
Least-squares extrapolation; Neural network 


